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Abstract. Recent differential measurements of elliptic flow are used to probe 
several hydrodynamic scaling predictions. Eccentricity scaling is observed for 
Cu+Cu and Au+Au collisions at ^snn = 200 GeV, suggesting essentially 
complete thermalization of the high energy density matter produced in these 
collisions. An estimate of the speed of sound is also obtained from the eccen- 
tricity scaled elliptic flow. The predicted mass scaling is observed for particles 
with transverse kinetic energy KEt up to ~ 1 GeV. For KEt values above 
~ 1 GeV, valence quark number scaling compatible with partonic degrees of 
freedom is observed. 
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1. Introduction 

Experimental measurements at the Relativistic Heavy Ion Collider (RHIC) have 
indicated the creation of matter having energy density in excess of what is required 
for the formation of a de-confined phase of quarks and gluons (QGP)[|T] |2 |3]. 
The decay of this matter results in particle emissions which show large anisotropics 
in their azimuthal distributions relative to the reaction plane. Such anisotropies, 
which are commonly characterized by the second Fourier coefficient v<i , suggest that 
the created high-energy density matter rapidly thermalizes and its decay dynamics 
is driven by substantial pressure gradients throughout its dynamical evolution [ 

HEDEJ. 

This picture for the matter created in RHIC collisions is supported by the 
good agreement achieved between i> 2 measurements [0 03 El El Ell E3 an( i those 
obtained form hydrodynamic model calculations [ 1131 1141 HB] . However, the detailed 
properties of the matter (ie. its equation of state, transport coefficients, effective 
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degrees of freedom, etc), are still lacking. Recently, several questions concerning 
the degree of thcrmalization of this matter have even been under debate [ 1151 116| . 
In this contribution, we show that the scaling characteristics of azimuthal 
anisotropy provide an important framework for developing important insights rele- 
vant to the properties of the high energy density matter created in ultrarelativistic 
heavy ion collisions at RHIC. 

2. Data & data analysis 

The results presented here utilize data obtained by the PHENIX collaboration for 
the Au+Au and Cu+Cu colliding systems at v / sjv r /v=200 GeV. A detailed descrip- 
tion of the PHENIX detector is given elsewhere [Ej- The procedures for tracking 
and event selection are described in Refs. [El EJ- The azimuthal anisotropies rel- 
evant to this work were measured using the reaction plane method [ 1181 1191 I2()j . 
This method measures the correlations between the azimuthal angle of the particles 
detected in the central arms of PHENIX and the azimuth of the reaction plane 
obtained with the two Beam-Beam Counters (BBC) located at | r\ |~ 3 — 3.9 [E]. 
An important advantage of the large r\ separation between the BBC's and the cen- 
tral arms of PHENIX is the expectation that the evaluation of the azimuth of the 
reaction plane is much less affected by possible non-flow contributions, especially 
those due to jets[E]- 

The magnitude of these elliptic flow correlations was determined by evaluating 
the second Fourier coefficient v<i , of the azimuthal distribution of emitted particles 
relative to the azimuth of the reaction plane [ 1181 1191 I2()| : 



l'2 



= ( e «(*-*«')) = (0062(0-*^)). (1) 



Here, <j> is the azimuthal angle of a particle measured in the laboratory coordinate 
system and §rp is an estimate of the azimuth of the reaction plane. A resolution 
factor was evaluated and used to correct for the dispersion of the measured reaction 
plane about the true reaction plane. 

Detailed measurements were performed over a broad range of centralities, pt 
and particle species. The results from these measurements are presented in the 
following sections. 

3. Insights into the degree of thermalization 

3.1. Eccentricity scaling 

As mentioned earlier, the large V2 values measured at RHIC [0 |H1 El Ell HH fTS] are 

compatible with those obtained from hydrodynamic model calculations [El 1141 ITfi| . 
Local thcrmalization constitutes an important assumption for the legitimate use of 
hydrodynamics to describe the data. This being so, a crucial question is whether or 
not one can gain further insights into the degree of thermalization from the data. 
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To do this, we recall here that an important prediction of the hydrodynamic model 
is that V2 should (i) scale with eccentricity and (ii) be independent of the colliding 
system for a given eccentricity [El El HH The eccentricity e = (y 2 ) — (x 2 )/((y 2 } + 
(x 2 )), measures the asymmetry of the "almond shaped" collision zone containing 
the high energy density matter created in the collision; x and y are, respectively, 
the short and long axes of the almond. 




Fig. 1. vijt as a function of c s (t — to)/R for different impact parameters. The plot 
is taken from [I15|. 



The eccentricity scaling prediction of hydrodynamics is illustrated in Fig. ^ 
The figure shows that, for a system of transverse size R, elliptic flow develops over 
a time scale ~ Rj (c s ) for matter with a speed of sound c s [ I15| . More importantly, it 
shows the variation of the scaled anisotropy V2/e, as a function of the characteristic 
distance c s (t — to)/R, where to is the thcrmalization time. The scaled results, which 
are shown for several impact parameters b or centralities, are taken from Ref. [I15j. 
They indicate that vijt is essentially independent of centrality except for the most 
peripheral collisions, where hydrodynamics is expected to break down. This result 
provides a clear signature for hydrodynamic scaling which can be explored with 
data. 

The left panels of Fig. show recently measured differential elliptic flow V2(pt) 
for several centralities in Au+Au and Cu+Cu collisions at ^/snn— 200 GeV. For 
both systems, the characteristic increase of V2(pt) as collisions become more pe- 
ripheral is evident. In a hydrodynamic scenario, this increase is expected because e 
grows as collisions become less central and larger values of e lead to larger pressure 
gradients which drive the magnitude of collective flow. The crucial question here is 
whether or not these data show the eccentricity scaling predicted by hydrodynamic 
models. 

The right panel of Fig.Elshows the scaled values of V2(pt) for the same central- 
ities shown in the left panel of the figure. To estimate the eccentricity, we use the 
fact that the p-r-integrated V2 values vary linearly with centrality and are essentially 
proportional to e, i.e e ~ k x V2 over a broad range of centrality selections. Here, 
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it should be emphasized that this proportionality has been experimentally verified 
in Au+Au collisions [0 |S|- From Glauber model estimates [[7] of e, we find that 
k ~ 3 for the cuts employed in this analysis. This factor is reflected in the right 
panel of Fig. [5] where we have used U2(p:r)/3w2(ceiitrality) for the scaled anisotropy. 
This ratio has the added advantage of being less affected by the systematic errors 
associated with our V2 measurements since they cancel out. 
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Fig. 2. The left panels show V2 vs. pt for several centrality selections for Au- 
(top) and Cu+Cu (bottom) collisions at ^/snn 



Au 



200 GeV. The right panel shows 
vi(jPt) scaled by the pr-integrated vi values for the same centralities indicated in 
the left panels. 



The scaled v^ values shown in the right panel of Fig. [21 indicate rather good 
scaling for the full range of centralities (or e) presented. This is clearly consistent 
with the eccentricity scaling predictions of hydrodynamics. Equally important is the 
fact that these results also confirm the expected independence on colliding system 
size. 

We interpret these findings as an indication that essentially full thermalization 
is achieved over a broad set of these Au+Au and Cu+Cu collisions [El \HZ\- The 
observation of the independence of v^/e on the system size also satisfies the require- 
ment of scale invariance which is a basic feature of perfect fluid hydrodynamics. 

Given these observations, it is important to ask whether or not one can use 
them to extract constraints for the properties of the high energy density matter. 
This is pursued in the next section. 



Scaling Characteristics of Azimuthal Anisotropy at RHIC 5 

3.2. An estimate of the speed of sound from scaled V2 data 

An important input for the hydrodynamic model is the equation of state (EOS). 
The EOS relates the pressure p to the energy density e. For a simple EOS, p = c s 2 e. 
Since v-i is driven by pressure gradients in the hydrodynamical model, the value of 
the speed of sound c s can greatly influence the magnitude of v? . 




Fig. 3. The left panel shows V2Je as a function of c s (t — to)/R for various values 
of c s . The calculation was performed for Au+Au collisions with impact parameter 
b = 8 £m at ^/sJTn — 200 GeV [El- The right panel gives a summary of the 
dependence of vije on c s . The dotted curve is represent a fit to the simulated data. 

The results from a recent study of the influence of c s on the magnitude of i>2/e [ 
1151 I2T] is shown in Fig. [31 The figure indicates that the magnitude of t^/e is rather 
sensitive to the assumed value for c s . This sensitivity is made more transparent in 
the right panel of Fig. |21 where we have summarized the results. 

An estimate of the speed of sound can be made by comparing the eccentricity 
scaled V2 values shown in Fig. |3 with the model predictions shown in Fig. |3 It 
is noteworthy that we have checked that the eccentricities used in the model are 
comparable to Glauber based model estimates [[Tj used in our evaluation of k. It 
should be pointed out that the v-ijt values in the left panel of Fig. differ from 
those in the right panel by a factor of 2. This factor was applied to take account of a 
well known difference in the procedure used to evaluate vi from the data and theory 
[I15j. A comparison between the right panels of Figs. [3 and gives the estimate 
c s ~ 0.35 ± 0.05 for a (px) ~ 0.45 GeV/c [ [22] ■ This value of c s is indicative of a 
soft equation of state [111]), with the matter possibly spending little time, during 
its evolution, in a region where c s = 0. 

3.3. Transverse kinetic energy scaling 

Given the observed eccentricity scaling, one can try to explore other hydrodynamic 
scaling scenarios. For instance, the variation of V2 with transverse kinetic energy 
KEy = rriT — rn, where itit = \/p\ + in 2 is the transverse mass of the particle, is a 



(i 



M.Issah and A. Taranenko for the PHENIX Collaboration 



good candidate for study. The use of this variable is intuitive because pressure gra- 
dients, which give rise to azimuthal anisotropy, lead to collective transverse kinetic 
energy of the emitted particles. 
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Fig. 4. Variation of V2 with (a) px and (b) KEt for selected hadrons for minimum 
bias collisions. The pion (it ), kaon (K ) and (anti-)proton ((p)p) results are from 
PHENIX [ I2H) and the neutral kaon Kq , lambda A and cascade 5 results are from 
STAR [EH. 



Figure 0] shows the dependence of V2 on pr and KEt (respectively) for several 
identified particles. The values of V2 for charged pions (7r ± ), kaons (K^) and (anti- 
)protons ((p)p) are preliminary data from the PHENIX Collaboration [ 1121 l2~3] 
whereas those for Kq, A and S are published data from STAR Collaboration [ 
I24| . The STAR V2 values were multiplied by a factor of 1.1 to account for a small 
difference in the mean centrality for minimum bias events measured by the two 
experiments. In the left panel of Fig. 0J one can observe mass ordering of i>2 (pt) 
for pt up to ~ 2 GeV/c. That is, for a given px, the larger the mass of a particle, 
the smaller its V2 ■ This pattern has been identified as a signature for hydrodynamic 
expansion [|2~5]. 

In the right panel of Fig.0] V2 is plotted as a function of KEt- Here, it can be 
seen that V2{KEt) shows the expected scaling with KEt up to ~ 1 GeV. For larger 
values of KEt, the scaling breaks. However, a rather striking trend emerges : the 
mesons and baryons scale independently for the whole KEt range. This pattern 
suggests that the constituents of which baryons and mesons are made also flow. We 
pursue this lead via an additional scaling test in the next section. 
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4. What are the degrees of freedom in the high energy density 
matter ? 

The observation in Fig.0]that mesons and baryons scale in two separate branches for 
KEt values above ~ 1 GeV suggests that the constituents of the hadrons flow. If this 
is indeed the case, then the data should reflect constituent quark number scaling. 
Indeed, several models based on the coalescence or recombination of constituent 
quarks have been developed [ 1261 1271 1281 1291 I30| . In these models, mesons and 
baryons obtain their elliptic flow from the elliptic flow of partons in an additive 
way, so that we have the simple relations: V2.m(pt) — 2 V2 iP (pt/2) and ^2,b(Pt) — 
3 W2, p (pt/3), where V2,m, V2,b and V2, P are the V2 of a meson, baryon and parton 
respectively. The underlying assumption of this picture is that each parton carries 
the same fraction of the hadron elliptic flow. 



0.1 



0.05 



(a) 



T 



T 



T 



• ji(PHENIX) Op (PHENIX) 
■ K (PHENIX) o A (STAR) 
K°(STAR) nE(STAR) 




I ' ' ' ' I 



(b) 



PHENIX Preliminary 




0.5 1 1.5 

p T /n q (GeV/c) 



0.5 1 1.5 

KE T /n a (GeV) 



Fig. 5. Variation of V2Jn q with pr/riq (left panel) and with KET/n g (right panel) 
for different hadron species for minimum bias collisions. 



The results from constituent quark number scaling are shown in Fig. \5\ Here 
the values for v%, pt and KEt for different particle species are divided by their 
respective number of quarks for mesons (n q =2) and baryons (n q =3). The left panel 
shows the commonly exploited quark number scaling ansatz: v%(jpT)/n q vs pT/n q . 
The left panel shows poor scaling for pr/n q < 1 GeV/c and better scaling, within 
errors, ior pr/n q > 1.3 GeV/c. In contrast, the right panel shows that V2(KEt)/?i 9 
scales over the whole range of KEx/n q . This observation can be interpreted as a 
strong indication that the pertinent degrees of freedom in the flowing high energy 
density matter are linked to the number of constituent quarks of the hadron. 
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Fig. 6. Variation of vijn q with KEr/n g for a comprehensive set of hadrons. 



Fig. shows the scaling with KEx/n q for a more comprehensive set of hadron 
species. The <f> meson and deuteron v-i values are taken from [|^ and [El respec- 
tively. It appears that the scaling holds even for a heavy meson like the <f> an d for 
the deuteron (n g =6), lending support to the idea that the different particle species 
emerge from a common flow held. 

4-1. Quark number scaling for different centrality selections 

Hydrodynamics is expected to break down for the most peripheral collisions and for 
relatively high pt[I13|. Consequently, it is important to investigate the evolution of 
the scaling picture as the centrality of collisions is varied and the px of the particles 
is increased. Such an analysis is still under intense investigation. Nonetheless, we 
show preliminary scaling results for a select set of centralities in Fig.0 Interestingly 
quark number scaling of v-i and KE-r holds rather well for the range of centralities 
presented. It will be interesting to see if this scaling breaks for more peripheral 
collisions and at high px ■ 

5. Conclusions 



In summary, we have presented results for the scaling properties of azimuthal 
anisotropy at RHIC. The eccentricity scaled v-i values for Cu+Cu and Au+Au at 
y'SNN = 200 GeV are similar as a function of px-i showing that scale invariance of 
azimuthal anisotropy is satisfied. This suggests essentially complete thermalization 
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Fig. 7. Variation of v%ln q with KE^/rig for 77*, if and p(p) for selected central- 
ities for Au+Au collisions at ^/snn = 200 GeV. 



of the high energy density matter produced in these collisions. The variation of v-i 
with transverse relativistic kinetic energy KEy shows the predicted hydrodynamic 
mass scaling at low KE^. However, it evolves into a differentiation of the mesonic 
and baryonic flow character of particles produced with high KEt- This suggests 
that KEt is a good variable to study the dynamics of the produced matter. Further 
scaling of v-i and KE^ by the number of valence quarks leads to a unified scaling 
behavior over the whole range of KEx/n q spanned by the data. The observations 
from eccentricity and transverse kinetic energy scaling point to the production of 
thermalized high energy density matter with partonic degrees of freedom. 
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